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Abstract 


An entropy correction method for the 
unsteady full potential equation is presented. 
The unsteady potential equation is modified to 
account for entropy jumps across shock waves. 
The conservative form of the modified equation 
is solved in generalized coordinates using an 
implicit, approximate factorization method. A 
flux-biasing differencing method, which 
generates the proper amounts of artificial 
viscosity in supersonic regions, is used to 
discretize the flow equations in space. 
Comparisons between the present method and 
solutions of the Euler equations and between the 
present method and experimental data are 
presented. The comparisons show that the 
present method more accurately models solutions 
of the Euler equations and experiment than does 
the isentropic potential formulation. 

Nomenclature 

speed of sound 
metrics of coordinate 
transformation. Equation (7) 
airfoil chord 
time step in computational 
space 

indices of grid points 
identity matrix 
Jacobian of coordinate 
transformation 
reduced frequency based on 
semi chord 
Mach number 
total speed 

distance from point of rotation 
to a point in the flow field 
gas constant 
entropy 
time 

contravariant velocities, 
Equation (6) 
cartesian coordinates 
angle of attack 

p 2 ' Y /J 

ratio of specific heats 
jump in potential across the 
wake 
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Subscripts 
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Superscripts 
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difference operator 

computational coordinate 

directions 

polar angle 

airfoil pitch angle 

density 

P/J 

biased density 

computational time 

velocity potential 

average of quantity across the 

wake 


isentropic 
normal to shock 
steady airfoil location 
pitch axis location 
points on airfoil 
minimum airfoil pitch angle 
maximum airfoil pitch angle 
free stream conditions 


lower side of wake 
iteration number 
upper side of wake 
sonic conditions 


Introduction 


Modern aircraft typically operate at high 
speeds where aeroelastic instabilities are more 
likely to occur. To successfully predict and 
analyze such aeroelastic phenomena, the aircraft 
designer needs methods that accurately predict 
the aerodynamic loads--steady and unsteady--that 
the vehicle experiences. Since many critical 
aeroelastic phenomena occur at transonic speeds, 
methods based on linear aerodynamic theory 
cannot accurately predict many aeroelastic 
responses. Thus, it Is necessary to use an 
aerodynamic method that can predict time- 
accurate solutions of nonlinear flows and that 
can accurately model shock waves and their 
unsteady motions. 

When shock waves appear in transonic flow 
fields, aerodynamic loads predicted using 
potential flow theory can be grossly inaccurate 
and even multivalued. Multiple solutions of the 
potential equation were first observed in two 
dimensions by Steinhoff and Jameson. 1 Salas 
and Gumbert2 showed that the phenomenon is not 
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confined to a particular airfoil or flow 
condition. Williams et al. 3 calculated 
multiple solutions using two-dimensional (2-0) 
transonic small disturbance (TSD) theory. 
Gibbons et al.4 showed multiple solutions for 
three-dimensional (3-D) TSO calculations on 
high-aspect-ratio wings. For lower aspect 
ratios, multiple solutions were not observed, 
but calculated lift coefficients were highly 
inaccurate when shock waves were in the flow 
field. 

Since potential theory can yield inaccurate 
transonic aerodynamic loads, aeroelastlc 
analysis performed using these loads has to be 
considered unreliable. The inaccuracy Is 
primarily a result of potential theory not 
modeling the jump In entropy that a fluid 
particle experiences as It passes through a 
shock wave. Therefore, calculated shock waves 
can have the wrong strength and be In the wrong 
location. 


p t + H> x + H) z = 0 (1) 

where $ is the velocity potential, and density 
P Is determined from Bernoulli's equation 

p = [1 + ^±(M* - 2* t - - ♦J ) f ?=T (2) 

The spatial coordinates, x and z, are normalized 
by airfoil chord c, and time t Is normalized by 
a„/c, where a„ is the free stream speed of 
sound. Density is normalized by the free stream 
density p«,, and $ is normalized by a«c. 

A transformation to a body-fitted 
coordinate system Is given by 

5 = ?(x,z,t) 

C = C(x,z,t) (3) 


Currently, the most widely used methods for 
unsteady, nonlinear flow analysis are based on 
TSD theory. 5-7 Fuglsang and Williams 8 
modeled the effects of entropy jumps through 
shock waves in 2-D TSD theory; Gibbons et al.4 
extended that method to three dimensions. These 
efforts resulted in TSD methods that more 
closely model solutions of the Euler equations 
than does isentropic TSD theory. 

However, TSD theory has some significant 
limitations. Only flows past bodies of small 
thicknesses at small angles of incidence and 
undergoing small amplitude unsteady motions can 
be analyzed. In addition, the TSD equations are 
formulated assuming that the free stream Mach 
number is near unity. Because of these 
limitations, the purpose of the present effort 
is to develop an improved method for predicting 
nonlinear, unsteady aerodynamic loads on 2-D 
bodies. A primary objective of this effort is 
to model the nonisentropic effects caused by 
shock waves. 

The present method is based on the unsteady 
full potential equation. Hafez and Lovell* 
presented a method for modeling entropy jumps 
across shock waves in steady potential flows. 
They showed that by modifying the isentropic 
density to include the effects of entropy jumps 
across shocks, steady potential solutions closer 
to Euler solutions are obtained. In the present 
effort, Hafez's and Lovell's method is extended 
to unsteady flows. 

A flux-biasing differencing methodlOjH 
Is used to discretize the flow equations and to 
model shock discontinuities. To increase 
computational efficiency, grids for unsteady 
calculations are generated, at each time step, 
using linear interpolation between grids 
calculated at extreme airfoil positions. 


x = t 

where 5 and c are the computational coordinate 
directions around and normal to the airfoil, 
respectively, and r is the computational time. 
Strong conservation form of (1) is maintained by 
writing the continuity equation in transformed 
coordinates as 


(4) 


If), * (f ) S ♦ &) t ■ 0 

Equation (2) transforms to 

1 

p = {1 +^i[Mf-2* T -(u+e t )* 5+ (w+c t )^]} r ' 1 (5) 

where the contravariant velocities in the £ and 
C directions, l) and W, respectively, are given 
by 


U = e t + A l 4> 5 + a 2 + c 

w = c t + a 2 * 5 ♦ a 3 * ? 

The metric terms are 

A, = 5 2 + e 2 

1 x s z 

A, = 5 C + £ r 

2 x x z z 


A = r 2 + c 2 
*3 ? x + h 


(6) 


(7) 


and the Jacobian of the transformation J is 
J - 5 X ? Z “ ^z ? x 
Approximate Factorization 


Problem Formulation 
Governing Equations 

The formulation used in the present effort 
is that presented by Bridgeman et al. 32 The 
flow field is described by the two-dimensional, 
unsteady full potential equation in conservation 
form 


Equation (4) is solved using first order 
backward differencing in time and second order 
central differencing in space as in reference 
12. The time derivative of density is 
linearized about previous time levels such that 
conservation form is maintained. The resulting 
equation is approximately factored into the form 
LgL ? di(i = F. This becomes 


2 
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Using (11) and (12), the flow tangency condition 
is imposed implicitly on both sides of (9) by 
changing the and L ? operators to 


x [I + hW n 6 c - ^r5 c (?A 3 ) n 6 c ](* n+1 - 4> n ) 


4 = [I + hU n 6 c - "(Aj - 


2 ,n , 

l - *7 } V 


(♦" . 4> n ~l J + A—(»" - 2d> n_1 + *"' 2 ) 

e 

+ h _( 5 " - p n -l) 
e n 


+ h ^n — (u n_1 <s c + w n_1 6 c )($ n - * n_1 ) 

+ ^[ 5 5 ( £ j) n + V 5W) " ■ F » ] 


L 5 = [i + hU n 6 5 - ^•S 5 (p A i) n « 5 J 


L? = [j + hW n 0 ? . ^{p-A 3 ) n 6 


= <j, n+1 . ,f, n 

and F represents the right side of (9). In (9), 
65 and 6 ? represent central difference 
operators, the superscripts n, n-1, and n+1 
represent time levels, 


■ l' 

and by representing the 6 t (pW) term on the 
right side as 2{pW)i .J+l/2* 

Dynamic Grids 

To apply the flow tangency condition at the 
instantaneous surface position of moving 
airfoils requires a new grid at each time step. 
Since the resources needed to compute these 
grids using an elliptic method 14 can be more 
than those necessary to do the aerodynamic 
calculations, an efficient interpolation method 
is used to generate the required grids. To 
simulate harmonic motions, the elliptic 
method 14 is used to calculate grids at the 
extremes of the motion. Grids for all other 
airfoil positions are then generated using 
Interpolation. Interior grid points are 
redistributed at each time step, while points on 
the outer boundary remain fixed. 

The method is illustrated in figure 1 for 
an airfoil pitching about a point Xp. A polar 
coordinate system centered at Xp is used. 
Using the subscripts 1 and 2 to denote the 
minimum and maximum pitch angles, the position 
of a grid point at any time x is given by 

9 s ( t) ' 9 sl 

r(x) - r. ^ T s — (r, - r.) (13) 

s2 si 


e (x) - e 

®( t )“ 8 i + T ' V (14) 

1 s2 9 sl 1 1 


The density biased in the £ direction is given 
by p, and F«. is a correction to the residual 
due to incomplete metric cancellation. > 1.3 

Airfoil Boundary Conditions 


where 0 s (x) is the instantaneous airfoil pitch 
angle, 


r l = (x l + z l )Z 


The flow is required to be tangent to the 
airfoil boundary. This condition is imposed by 
requi ring 

W = C t + Ag* e + A 3 $ c = 0 (10) 

at the airfoil. Equation (10) is applied as 

<*C>1.J " ' ( ‘S7 *e } i,J ' ^i.J (U) 


where the subscript J represents points on the 
body. The tangency condition implies that 


- (4 + 4)‘ 


e = tan _1 (— i) 
x Xj 


9 2 " 


^1, 0-1/2 = " ^1, J+l/2 


Interpolated grid points are then given by 
x(x) = r(x)cose(x) 
z(x) = r(t)sin0(x) 
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Points on the body, x s and z s , are defined 
at each instant by 

x s ( t > " x p + (V x p) cos 6 s( T ) + z o s 1 n V T > 


Flux Biasing 

The residual terms on the right side of (9) 
are centrally differenced about the node point 
(1,j) to give 


z s (t) = - (x Q - x p )s1ne s {t) + z 0 cose s (r) 

where the subscript o represents steady state 
locations. The time-dependent metric terms 


are then computed using first order differences 
for x T and z x . 

To demonstrate the interpolation method, 
figure 2 shows grids for an airfoil pitching 
about its quarter chord and the angle of attack 
<*( t) given by a(x) = 22.5°(1 + sin(k t)). The 
calculated grids for a(x) = 0° and 45° are shown 
in figures 2a and 2d, respectively; they were 
generated with the GRAPE program 1 ^ which uses 
an elliptic method. The interpolated grids for 
a(x) = 15° and 30° are shown in figures 2b and 
2c. Even for such large angles, the 
interpolation works successfully. The details 

of the near-field grids cannot be seen, but the 
grids remain near orthogonal; the metric A2, 
which is a measure of the skewness of the grids 
is much smaller than Aj and A3 over the 
entire grid. If A2 = 0 everywhere, the grid 
would be perfectly orthogonal. Grid lines 

intentionally were not plotted across the wakes 
to make those regions easy to identify. 


S „ (£Ii\ t PU\ 

VJ ; i,j 'J M+1/2.J '1-1/2, j 


. /PU 


pIK 


M pW) i,j = (pW) i,j+l/2 ' (pW) i , j-1/2 


The computational cell is shown in figure 3. In 
supersonic regions, artificial viscosity, 
necessary to capture shocks, is introduced by 
biasing the density in' the upwind direction. 
Although it is generally necessary to upwind the 
density in both computational directions, in 
this work, upwinding only in the £ direction has 
been fo^nd to be satisfactory. The biased 
density p is defined as (for U > 0) 


1+1/2, j q i+ i/ 2 j 


-[pq - A5(pq)'] 1+1/2J 


q 1+l/2,/ (pq)i+1/2 ’ j " (pq)l+ l/2,j 


where 

+ {pq) i-l/2,jJ 

(16) 

(pq) - 

= pq - p*q* q > q* 

0 q < q* 

(17) 


Unsteady Wake Condition 

For lifting flows, the shed vorticity is 
represented as a jump in potential across a wake 
line. The following wake boundary condition 12 
is used in the present effort 

r x + <w>r = 0 (15) 

where r is the jump in potential across the 
wake, <(> u - and <W> is the average of VI 

above and below the wake, 1/2 (W u + W*). The 
unsteady transport condition for the jump in 
potential, (15), is obtained by assuming that 
the isentropic density and normal velocity are 
continuous across the wake. It does not model 
the jump in entropy across the wake that can 
result when shocks of different strengths are on 
the upper and lower surfaces. However, Hafez 
and Lovell q have demonstrated that assuming 
the isentropic density to be continuous across 
the wake is a good approximation for the full 
potential formulation with entropy corrections. 

Far-Field Conditions 

In the far field, the flow is set to free 
stream conditions 

♦ = H x 

' 00 


poo 

q is the flow speed (q = + $+), and q* and 

P* are the sonic speed and density, 
respectively, given by 

+ I r( M -- 2 *x- 2 «tV 2 <t* c )J 

(IS) 

P* = (q* 2 )^ (19) 

For steady flows, p* and q* are constants that 
are computed once, but for unsteady flows, p* 
and q* must be computed at each grid point at 
every time step. 

Flux biasing is an improvement upon density 
biasing since it (a) accurately tracks sonic 
conditions and requires no empirical constants 
to specify the amount of artificial viscosity, 
(b) produces no velocity overshoots at shock 
waves, thus allowing for larger time 
steps — increasing computational eff 1ciency--for 
unsteady calculations, (c) produces well 
defined, monotone shock profiles with a maximum 
two point transition between the upstream and 
downstream states, and (d) dissipates expansion 
shocks, ruling out solutions with such 
nonphysical characteristics. Shankar!5,16 has 
previously used flux biasing in unsteady full 
potential calculations. 


p = 1 

where H» is the free stream Mach number 
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Entropy Correction Method 


At shock points, (23) is modified to 


When a fluid particle passes through a 
shock wave, it experiences a jump in entropy. 
The entropy jump as is a function of the normal 
Mach number upstream of the shock M n 


AS 

R 


1 ilnf 2Y M 2 yin [ (Y+1)Mn 

FT tln( FT M n " FP ' y1 " L 


(y-i)m; + 2 


-]} 


(20) 


-AS 


-AS 


r #- _ R \n+l ^ R xn il 

L( p 1 e )i,j - < p i e >i,j>F 


n+1 

fPM.i 


n+l 


+ . fPLLi 

'1 +1/2, j 'J '1 -1/2, j 

-AS -AS 

+ (pe R w )i , j+l/2 ' (pe R W) i,j-l/2 = 0 


( 24 ) 


For unsteady flows, the shock speed must be 
monitored to obtain the speed of the flow 
relative to the shock. Here, the entire shock 
is assumed to move at the same speed as the base 
of the shock. To account for the entropy jump, 
the density is modified to 


where as is a function of M ni 3 ^ 2 y and 
-As 

p i +1/2 , j = (p i e J i+l/2,j + 


- AS 

P = P^ R (21) 

where the isentropic density, pj, is given by 
(5). Expanding the continuity equation, (1), 
yields 


-AS 

e R h t + ( p i*x>x + (p i*zU - 


a ( p i0 - p *0*)f . 1 /? i 

q i+l/2,j 1 1 1/Z,J 


(25) 


The computational cell at shock points is shown 
in figure 4. In the approximate factorization 
scheme, this is implemented by representing 


the iL(p n - p n_1 ) and ^[6^) + « (pH)] 
terms in (9) as 


»1* 'Hf-lt (f)zl-0 


( 22 ) 


-AS 

R 


AS 

R 


V (p i e )i,j-( p i e >?;]]•"< 


where u = $ x , and w = <f> z . The last part of 
(22) vanishes since ^ = o in the isentropic 

flow regions. The remaining part of (22) is the 
classical potential equation except across shock 
waves and wakes where there is a jump in 
entropy. Hafez and Lovell q showed that 
assuming the isentropic pressure to be 
continuous across the wake is a good 
approximation except near the trailing edge. 
Thus, that assumption is made here, and the wake 
condition in (15) is used for all nonisentropic 
calculations. 


II 11 II 11 

n | / pU x /pUx 

B n u J ^ i +1/2, j ’ M-l/2, j 


-AS -AS 

+ (p i e R W) i, j+l/2 ‘ (p i e R W) i,j-l/2^ 


respectively, where Pj +1/2 , j is given by (25). 
The effect of the entropy correction is to 
introduce a source distribution along the 
shock. The source strength is dependent upon 
AS. 


To implement the entropy corrections, the 
difference operator at shock points is 
modified. Equation (4) is discretized as 


n+l - n % 

'f.fi ■ 1 » j + (£“ 

d 1+1/2, j 



n+l 

i -1/2 ,j 


+ (pW) 


n+l 

1, j+l/2 


(pH) 


n+l 

i.j-1/2 


= 0 


(23) 


Results and Discussion 


Steady pressures on an NACA 0012 airfoil 
have been calculated using the isentropic and 
nonisentropic formulations. Euler 

calculations 17 at the same flow conditions 
were also made. Unsteady pressures on the NACA 
0012 oscillating in pitch about its quarter 
chord were computed using the isentropic and the 
nonisentropic methods. Comparisons of those 
calculated pressures with TSD calculations and 
with experimental data were made. The above 
unsteady calculations were made using grid 
interpolation to simulate the airfoil unsteady 
motion. In addition, unsteady calculations were 
made using grid rotation to simulate the airfoil 
motion. Comparisons of the unsteady results 
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obtained by the two methods were made. Steady 
pressures on the NLR 7301 airfoil were computed 
using the isentropic and the nonisentropic 
methods. Comparisons between the Isentropic 
method, TSD, and experiment and between both 
full potential methods, Euler calculations, and 
experiment were made. All calculations were 
made using 0-type grids with 101 points around 
the airfoil and 31 points outward from the 
airfoil. 

NACA 0012 

In this section, the calculations made for 
the NACA 0012 airfoil are presented. These 
results demonstrate the accuracy of the present 
Isentropic potential method and the improvements 
that result when the effects of entropy are 
modeled. 

Figure 5 shows the steady pressures that 
are calculated using the potential methods and 
an Euler method*? for = 0.84, a = 0°. 
For this airfoil, the flow conditions are in the 
region where multiple solutions are known to 
occur. The Isentropic result in figure 5 shows 
an asymmetric solution with negative lift. The 
other possible solutions are an asyirenetric 
solution with positive lift (the opposite of the 
solution in figure 5), and a symmetric solution 
with zero lift. The Euler solution for this 
case is symmetric with zero lift. When 

entropy corrections are used, the potential 
method yields a symmetric, nonlifting solution. 
The calculated pressure distribution agrees very 
well with that obtained by solving the Euler 
equations.*? By modeling the shock-generated 
entropy, the nonphysical asymmetric solution is 
eliminated, and a symmetric pressure 
distribution is obtained. 

Calculated unsteady pressures on an 
NACA 0012 oscillating in pitch about its quarter 
chord, a( t) = 0.016° + 2.51°sin(kT), at 

M«, = 0.755 and k = 0.0814 are shown in figures 
6a-6f. Included in the figures are isentropic 
full potential (FP) calculations, TSO 
calculations made using the code of reference 6, 
and experimental data.* 8 The instantaneous 
pressure distributions show that the present 
method agrees very well with experiment over the 
entire range of unsteady motion. These 

calculations show the Increased accuracy that is 
obtained by going from TSD theory to full 

potential and applying the flow tangency 
condition on the actual airfoil boundary. The 

most notable improvement is the accuracy of the 
computed pressures on the forward portion of the 
airfoil. The present method does a much better 
job of calculating the flow over this portion of 
the airfoil than does the TSD method. The FP 
shock positions and strengths show good 
agreement with the experimental data, although 
the shock waves generally are further aft and 
stronger than the measured shocks. This is the 
expected result for an isentropic potential flow 
method. 

For improved accuracy of the FP method, the 
effects of shock-generated entropy are included 
in the calculations. The corresponding results 
for the case of figure 6 are shown in in figures 
7a-7f. Those figures show the isentropic and 
nonisentropic pressures and experimental data. 
Generally, at each instant in the cycle of 


motion, the effects of the entropy corrections 
are to weaken the shock and move it forward such 
that the nonisentropic pressures agree better 
with the measured data than do the isentropic 
pressure distributions. In figure 7d, the 
position and strength of the isentropic shock on 
the lower surface are in excellent agreement 
with the measured data, and the effects of 
entropy cause virtually no change. Including 
the effects of entropy generally does not result 
in more accurate modeling of the flow 
immediately downstream of shock waves. This is 
because the thickening of the boundary layer 
that occurs downstream of embedded shocks is not 
modeled in the inviscid formulation. At points 
in the cycle where the shocks become strong, the 
measured pressures immediately behind the shocks 
show the effects of boundary layer thickening. 

In calculating the pressures of figures 6 
and 7, the airfoil motion was simulated by 
Interpolating the grids between the extreme 
airfoil positions. An attempt was made to 
simplify the calculations by rotating the entire 
grid to model the airfoil motion. The metrics 
were recomputed at each time step, and the 
airfoil boundary condition was applied at the 
instantaneous surface position. Examples of the 
unsteady pressures calculated using this 
procedure are shown in figures 8a and 8b. The 
calculations correspond to those in figures 6a 
and 6e, respectively. Included in the figures 
are experimental data and the isentropic 
pressures calculated by interpolating grids to 
simulate airfoil motion. It is seen that 
rotating the entire grid causes large 
differences in the calculated pressures and very 
poor agreement with the measured data. Thus for 
the current formulation, grid rotation should 
not be used. It should be noted that Malone and 
Sankar* 9 successfully used grid rotation and 
translation for unsteady full potential 
calculations. That method uses a different 
linearization and different far-field boundary 
conditions than those used here. 

NLR 7301 


Steady pressures on the NLR 7301 airfoil at 
its design condition, M*. = 0.721 and 
ct = -0.19°, are shown in figure 9. Included are 
isentropic FP calculations, TSD calculations, 
and measured pressures. 2° Again, the improved 
accuracy that is obtained with the full 
potential formulation is evident. Although the 
level of pressures on the upper surface is 
higher than the experimental values and the 
shock is too strong and aft of its experimental 
location, the present method correctly predicts 
the trend of the pressure distributions. The 
TSD method predicts two shocks on the upper 
surface— one on the forward portion of the 
airfoil and another near the measured shock 
location. 


Entropy corrections were applied to the 
case of figure 9; the results are shown in 
figure 10. That figure shows a comparison of 
the isentropic and nonisentropic FP pressures, 
Euler*? pressures, and experimental data. 20 
The effects of modeling the shock-generated 
entropy are to weaken the shock and move it 
forward. The level of pressures upstream of the 
FP shock remain higher than those of the Euler 
pressures and the experimental data. 


6 



Concluding Remarks 

A method for modeling the effects of 
shock-generated entropy has been developed for 
the unsteady full potential equation. This was 
accomplished by modifying the isentropic density 
immediately downstream of shock waves to account 
for the local jump in entropy. The entropy 
correction method was implemented in a 
two-dimensional full potential code that used a 
flux biasing differencing method to discretize 
the flow equations. The resulting method was 
tested for steady and unsteady flow past an 
NACA 0012 airfoil and for steady flow past an 
NLR 7301 airfoil. Comparisons of the isentropic 
and nonlsentropic calculations were made with 
transonic small disturbance and Euler 
calculations and with experimental data. 

For the NACA 0012 airfoil, modeling the entropy 
effects alleviated the phenomenon of multiple 
solutions in a steady flow case. In that case, 
a unique solution that showed good agreement 
with an Euler solution was obtained. For 
unsteady flow, the primary effects were to 
weaken the shock and to cause slight changes in 
its location. In comparisons of the unsteady 
calculations with experimental data, the 
nonisentropic shock always moved toward the 
measured shock location. 

Calculations for the NLR 7301 airfoil at 
its design condition also showed the effects of 
the entropy corrections to be a slight weakening 
and forward shift of the shock from its 
isentropic location. The calculated pressures 
show the correct trends, but the pressure levels 
are too high. 

To summarize, the entropy correction method 
alleviated the phenomenon of multiple solutions 
in a case where it was observed. The primary 
effects of modeling shock-generated entropy are 
to weaken and cause a slight shift of the shock 
position. 
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Figure 2. Example of grid Interpolation. 
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a) cc(t) = 2.01°, kx = 127° 


d) q(t) = -2.41°, kT = 255° 



O Experiment upper 



b) c.(t) = 0.52°, kT = 168° 


e) cx(t) = -2.00°, kT = 307° 


O Experiment upper 



x/c 


O Experiment upper 



c) cx(t) = -1.25°, kT = 210° 


f) ct( t) = -0.54°, kT = 347° 


Figure 6. Isentropic unsteady pressure vs. experiment for an NACA 0012 airfoil 
oscillating in pitch, M. = 0.755, o(t) = 0.016° + 2.51°sin(kT). 
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a) ci(t) = 2.01°, kr = 127° 


d) o(t) = -2.41°, kT = 255° 



x/c 


b) o(t) = 0.52°, kt = 168° 



x/c 


e) o(t) = -2.00°, kx = 307 c 



x/c 


c) cx(x) = -1.25°, kt = 210° 



x/c 


f) o(t) = -0.54°, kx = 347° 


Figure 7. Isentropic and nonisentropic unsteady pressure vs. experiment for an 

NACA 0012 airfoil oscillating in pitch, M» = 0.755, q(t) = 0.016° + 2.51°sin(kx) 
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O Experiment lower 



x/c 



a) a(x) « 2.01°, kx = 127 0 


Figure 9. Isentroplc steady pressures vs. 

experiment for an NLR 7301 airfoil, 

M = 0.721, a = -0.19°. 
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O Experiment lower 
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Experiment upper 
Experiment tower 
Isentroplc FP upper 
Isentroplc FP lower 
Nonisentropic FP upper 
Nonisentropic FP lower 
Euler upper 
Euler lower 


b) a(x) • -2.00°, kx = 307° 


Figure 8. Unsteady pressures on an NACA 0012 

airfoil calculated using interpolated 
and rotated grids, = 0.755, 

a(x) = 0.016° + 2.51°sin(kx) . 


Figure 10. Isentroplc and nonisentropic steady 
pressures vs. experiment for an 
NLR 7301 airfoil, * 0.721, 

a = -0.19°. 
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